Indifference Pricing of American Option Underlying Illiquid 
Stock under Exponential Forward Performance 

Xiaoshan Chen * Qingshuo Song < Fahuai Yi * George Yin § 



Abstract 

This work focuses on the indifference pricing of American call option underlying a 
non-traded stock, which may be partially hedgeable by another traded stock. Under 
the exponential forward measure, the indifference price is formulated as a stochastic 
singular control problem. The value function is characterized as the unique solution 
of a partial differential equation in a Sobolev space. Together with some regularities 
and estimates of the value function, the existence of the optimal strategy is also ob- 
tained. The applications of the characterization result includes a derivation of a dual 
representation and the indifference pricing on employee stock option. As a byproduct, 
a generalized Ito's formula is obtained for functions in a Sobolev space. 
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1 Introduction 



Recent developments of indifference pricing, particularly the valuation of American options 
under forward performance measures, have attracted much attention; see [10] and the ref- 
erences therein. One of many applications of the indifference pricing is the valuation of 
Employee Stock Option (ESO); see [12, 11]. An ESO is an American type of call option on 
the common stock of a company, which is not available for trade in the market. However, 
there may be some other partially correlated stocks (ex. Index) available in the market, 
which may be used for the purpose of partial hedge of ESO. Such a problem can be formu- 
lated as a stochastic singular control problem [14] with optimal stopping. Related study 
also appears in empirical research [2, 7] for the study of early exercise behavior. 

In this paper, we study the indifference pricing of the American call options underlying 
a non-traded common stock by assuming that there exists another perfectly liquid stock 
available in the market, which may be used for the purpose of partial hedge. The indiffer- 
ence pricing involved is associated to a portfolio optimization problem, which is in fact a 
stochastic control problem with optimal stopping involved. Using the classical verification 
theorem [18], the value function may be characterized as the unique solution of an associ- 
ated variational inequality (VI) under rather strong assumptions. These assumptions are: 
(1) The VI has the unique classical solution, i.e., unique solvability in C 2 ' 2,1 ; (2) The control 
problem has an underlying process and associated control process that attains the optimal 
value, i.e., existence of the optimal pair (tt*,t*). We refer to [4, 18] for the verification 
theorem in general control theory, and [10] for the case of indifference pricing of American 
options. 

In the general setting of control problems, the aforementioned conditions for the veri- 
fication theorem are difficult to verify due to the full non-linearity. As an alternative, one 
can use the viscosity solution approach [3] to identify the value function. In particular, [15] 
characterizes the value function V of an investment problem in the framework of indifference 
pricing. A drawback of this approach is that due to the lack of information on the regularity 
for the viscosity solution, one cannot obtain further knowledge of optimal controls. Even 
the existence of the optimal control is problematic. 

In this work, we will adopt an intermediate methodology between the classical verifica- 
tion theorem and the viscosity solution approach. First, we show the verification theorem 
holds under weaker assumptions than the classical verification theorem needs, but stronger 
than the viscosity solution needs. For the purpose of complete characterization, we obtain 
some regularity estimates of the partial differential equation (PDE), which verifies all the 
necessary assumptions of the existence of the PDE solution and associated optimal pair. 
Thanks to the above results on control problem, one can further develop more useful re- 
sults: They are the regularity of optimal exercise boundary, the dual representation, and 
monotonicity of the value function on the system parameters. To be more precise, the idea 
is illustrated below. 

In the standard arguments of the verification theorem, one uses Ito formula and classical 
solutions of a variation inequality (VI). However, it is too restrictive, since the VI in our case 
does not admit a classical solution. Therefore, we replace the assumption on the existence 
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of classical solution by solutions in Sobolev space W p ' l( ^ c for large enough p. Then, the 
question is 
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• How large p is enough to enable us to use ltd formula? 

We will demonstrate that p = 3 is large enough. As a by product, we generalize Ito's formula 
given in [8] from one-side inequality to two-sided equality; see the proof in the appendix, 
which may be potentially useful in other related control problems. 

Moreover, utilizing special structure with the exponential forward performance, one can 
derive the verification for the price P associated to a PDE. Note that another assumption 
needed for the verification theorem, namely the existence of optimal pair, can be reduced 
to an estimation on the first derivative. As a result, the verification theorem holds under 
weaker assumptions on solvability (HI) and estimates (H2) given at the end of Section 3. 
To fully characterize the control problem, we answer the questions: 

• Does solvability (HI) and estimation (H2) hold for the VI associated to the price P? 

To proceed, we start with a simple transformation on the backward equation (in time), which 
leads to an equivalent counterpart of forward equation on standard domain Q = Rx (0, T]. 
Since the original domain is unbounded, we shall use penalized method on the truncated 
version. As a result, the truncated PDE after penalization is quasi-linear, and Leray- 
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Schauder fixed point theorem provides the solvability in W p ' , and the comparison principle 
is standard for the strong solution, see more details in Section 4. By forcing the limit 
of the parameter of the penalized function e — > 0, and the parameter of the truncated 
domain N —?■ oo, we use local estimate to obtain the C a,a l 2 estimate (De Giorgi-Nash- 
Moser estimate) and Wp' 1 estimate on compact domain by removing out singularity point 

2 1 

(In if, 0). Consequently, we obtain the existence of W ' loc (Q) nC(Q) solution with some 
regularity estimates on the first order derivatives, which eventually resolves (H1)-(H2). To 
this end, we complete the characterization, which is summarized in Theorem 11. 

One application of the main result is that, a dual representation Proposition 12 is proved 
in the framework of stochastic control theory. A simple derivation from this representation 
implies several economically interesting facts: sensitivity properties with respect to risk 
aversion 7 and some other parameters, and sublinear property with respect to the payoff, 
etc. Another application is the indifference pricing of ESOs, where vesting periods and 
job termination risk is involved. Thanks to the characterization theorem, the price can be 
characterized as the classical solution of its associated PDE. 

The rest of this paper is outlined as follows. The precise formulation of the problem is 
given in Section 2. Section 3 proceeds with a generalized verification theorem. Section 4 
presents the main results that gives a complete characterization. Section 5 is concerned 
with the dual representation and other properties. Section 6 demonstrates application to 
ESO costs. Finally, an appendix containing a few technical results is provided. 

2 Problem Formulation 

Let (f2, J 7 , P) be a filtered probability space with a filtration (J-t)t>o that satisfies the usual 
conditions, and the processes B and B be two independent standard Brownian motions. 
Let 7t,T be the set of all stopping times in the interval [t,T]. 
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We consider a financial market with a risk-free interest rate, which consists of a non- 
traded stock issued by a firm with its price Y„ v at time v € [t,T] given by 



dY v = bY u du + cY u {pdB v + y/1 - P 2 dB v ), Y t = y, (1) 

and one liquidly traded stock with its discounted price St' s at time v 6 [t, T] 

dS v = S u a(Xdv + dB u ), St = s. (2) 

The payoff of American call option with strike K > and maturity T underlying non- 
traded asset Y v ' v is g(Y T ), where r is the exercise time r £ 7t,T and g is given by a function 
given by 

g{y) = {y-K) + . 

If an employee, with initial capital w, dynamically trades in the stock S, then her wealth 
process W t ' w ' 7T under self-financing rule satisfying 

dWl = Ti v o(\dv + dB v ), W{t) = w, (3) 

t S 

where ir u represents the cash amount invested in the liquid stock SV . We assume the 
strategy it belongs to Zt, which is the set of all progressively measurable processes it : 
[t, T] x Q — >• R satisfying integr ability condition 



E 



T 



TT^df 



< oo. (4) 



We introduce an exponential forward performance measure given by 

U t {w) = - e -r»+l> 2 *; (5) 

see [10, Theorem 3 and (18)]. Under the above performance measure, we compare the 
following two scenarios for an employee holding the initial capital w and a unit of American 
call at the initial time t: 

1. Find a stopping time r 6 % t to exercise the call option, while dynamically trading the 
capital w in liquid stock S v , to maximize its performance, and find the corresponding 
value 

V(w,y,t)= sup E[U T (W^ + g(Y^))\T t ). (6) 

2. Receive a cash payment of amount P(w,y,t) by selling one unit of call at time t, the 
performance corresponding to the total capital w + P(w, y, t) is Ut{w + P(w, y, t)). 

Now, we are ready to define the indifference price of the American call option. The 
indifference price of this call is the cash payment P(w,y,t), which makes the above two 
scenarios indifferent with respect to the given forward performance of (5). Therefore, the 
price P(w,y,t) is the value satisfying 

V(w,y,t) = U t (w + P(w,y,t)). (7) 

To determine the above indifference price of the American call, we consider the following 
two problems: 
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(Ql) Find the indifference price using the identity P(w,y,t) = U t (V(w,y,t)) — w, after 
solving the stochastic control optimization (6), i.e., 

(a) characterize the value function V of (6); 

(b) characterize the pair of optimal strategy r* and ir* if it exists. 

(Q2) Characterize directly the indifference price P(w, y, t) as the unique solution of certain 
PDE. 

To proceed, the following assumptions are imposed throughout this paper: 
(Al) Assume a, 7 > 0, and p G [0, 1). 

We exclude the case p = 1, where the market is complete and the problem of indifference 
price can be reduced to a Black-Scholes model. 

3 Generalized Verification Theorem 

We present a version of the verification theorem for the value function V of (6) and P of (7) 
in this subsection, which generalize the verification theorem from the classical PDE solution 
to the strong solution in Sobolev function space with certain regularity. 

Define a three dimensional domain by Q = Kxl+x [0, T), and denote by Wffa '(Q 1 ), 
the collection of all functions Q} 3 (w,y,t) i-> ip(w,y,t) having (d ww (p,d yy ip,dt(p) in dis- 
tribution sense and integrable in any compact subset of Q . For a scalar ir € M, define a 
parameterized differential operator L\ : W^f^. ^ R as 

CJf(w, y, t) = ^c 2 y 2 d yy (p + byd y (p + ^a 2 7r 2 d ww f + ir(pacyd wy cp + \ad w ip). 

We also need to introduce the continuation region by 

C[v] = {(w, y, t) G Q 1 : v(w, y, t) > U t (w + g(y))}. (8) 

Lemma 1 (Verification Theorem for Function V) Suppose there exists v € W 2 f^ c (Q l ) 
for some p >3, satisfying 

min{ - d t v + inf {-C\v (w, y, t)}, v(w, y, t) - U t (w + g(y))} = 0, (9) 
with the terminal condition 

v(w,y,T) = U T (w + g(y)). (10) 

Then, v(w,y,t) > V(w,y, t). If, in addition, there exists a pair (W*,n*) satisfying (3), (4), 
and 

(d t v + £pv)(W;,Y v , v) = 0,Vt<v<T*, (11) 

where 

t* :=mf{v>t: (W*,Y u ,u) (£C[v]}AT. (12) 

Then the variational inequality (9)-(10) admits a unique solution in W^f^^O 1 ), and v(w, y, t) = 
V(w,y,t). 
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PROOF: Fix an arbitrary stopping time r € 7~t,T and an admissible control tt S Ztn an d 
denote (W,Y) := (W t ' w ^, ) and E t [ ■ } := E[ • \T t ] for simplicity. By virtue'of the 
generalized ito's formula (see Proposition 15), we have 

E t [v(W T ,Y T ,T)\ = v(w,y,t) +E t [£ (d t v + C\v)(W v ,Y v ,v)dv . 

Applying (9), we have v(w,y,t) > Ut(w + g(y)) for all (w,y,t), hence 

E t [v(W T ,Y T ,T)} >E t [U T (W T + g(Y T ))}. 
Note that (9) also implies (dtv + Civ)(w,y,t) < for all (w,y,t), which yields 



E t 

Therefore, we conclude 



(dtv + £[v)(W v ,Y v ,v)dv 



< 0. 



E t [U T (W T + g(Y T ))]<v(w,y,t), 
and arbitrariness of r and tt further implies one-sided inequality 

v(w,y,t) > V(w,y,t). 

Furthermore, if (11) holds for some (W*,tt*), Ito's formula together with the definition of 
C[v] gives the inequality of the opposite direction: 

v(w, y, t) = E t [v(W T * , Y r * , t*)] = Et [U r * (W T * + g(Y T *))] < V(w, y, t). 



This implies v = V. □ 

Next, we discuss the price P of (7). To proceed, we introduce a domain Q 2 

7 2 f ( 



[0, T) and its Sobolev space 

w vL(Q 2 )- Also defme 

a non-linear differential operator C2 ■ 



<L(2 2 )^^by 

£2<f(y, t) = ^c 2 y 2 dyy(p + (by - pc\y)d y (p - -7(1 - p 2 )c 2 y 2 (d y ip) 2 . 

Lemma 2 (Verification Theorem for Function P) Suppose there exists a function f £ 
Wpf oc (Q 2 ) for some p > 3, satisfying 

mm{-d t f(y,t)-C 2 f(y,t),f(y,t)-g(y)} = 0, (y,t)eQ 2 , (13) 
urai/i t/ie terminal condition 

f(y,T) = g(y), y G ]R + , (14) 

and |5 y /(y,i)| is uniformly bounded. Then, PDE (9) together with the terminal condition 
(10) is uniquely solvable in W f^{Q l )i an d there exists a pair (W*,tt*) satisfying (3), (4), 
and (11). As a result, the price function P(w,y,t) of (7) is independent to the initial wealth 
w, and P(w,y,t) = f(y,t). 
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Proof: Let v(w,y,t) = U t (w + / (y , t)) . Then, we have v € W^Q 1 ). One can directly 
check d ww v = j 2 v < in Q 1 , and hence, the function ir*[v] : Q 1 i— >• M of the form 



K*[v](w,y,t) 



pcyd wy v(w, y, t) + Xd w v(w, y, t) 



ad ww v(w,y,t) 

is well defined. Note that, ir*[v] is independent to the variable w, i.e., one can rewrite 

7r *[v](y,t) = -^.yd y f(y,t) + —. 

a <77 

Now, we can define a process it* by 

K = ^[v](Y u ,u) = -^-Y u d v f(Y u ,u) + —, u€(t,T). (15) 
a (77 

Note that, F is the unique strong solution of (1). Therefore, ir* is an admissible strategy, 
since it satisfies the integrability condition (4) due to boundedness of d y f, i.e., 









- r T 




E 




< C + CE 


it Y " dv 


< oo 













Since the optimal strategy ir* is independent to W*, the solution to (3) can be simply given 
by its integral form of (3) . Calculations by change of variables leads to (9)-(ll). Thus, 
Lemma 1 together with definition (7) implies V(w,y,t) = Ut(w + f(y,t)), and f(y,t) = 
P(w,y,t). □ 

Remark 3 (The boundary condition of f(y,t) at y = 0) The verification theorem can 
be thought of the probability counterpart of uniqueness result of PDE solution (13)- (14). 
From the above verification theorem, the uniqueness holds without the specification of the 
boundary condition on y = 0. Similar observations are addressed in [16] (Fichera condition) 
for linear equation, and in [1] for the fully non-linear parabolic equation without obstacle. 

To this end, we summarize the implication of Lemma 1 and Lemma 2. The above 
verification theorems give conditional characterization of the control problem since they 
give the answers to (Q1)-(Q2) based on the assumptions on the solvability of PDE. More 
precisely, 

1. For (Ql-a), V is the unique solution of PDE (9)-(10); 

2. For (Ql-b), the pair of optimal control exists, and they may have representation of 
(15) and (12), respectively; 

2 1 

3. For (Q2), P is invariant in variable w, and the unique W p ' loc solution of PDE (13)-(14) 

function of two variables of (y,t), 

provided that the following hypothesis are valid, 
(HI) (13)-(14) is solvable in W^ 0C (Q 2 ); 
(H2) \d y f(y,t)\ is uniformly bounded. 
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4 Main Result: Complete Characterization 



To obtain the complete characterization, it is crucial to study the solability and its related 
estimations of PDE (13)-(14). 

For the convenience in the analysis, we analyze the backward equation (in time) (13)- 
(14), by studying its associated forward equation of the following form. 

Let x = Iny, 8 = T — t, u(x, 9) = f(y, t) in (13)-(14), then u(x, 9) satisfies 



J mm{(d e u- Cu)(x,9), u - (e x - K)+} = 0, (x, 9) G Q, 
\ u(x,0) = (e x -K)+, i£K, 
where Q = R x (0, T] and £ is the differential operator given by 

C(p(x,t) = \c 2 d xx ip + (b - pcX - \c 2 \d x tp - -7(1 - p 2 )c 2 (d x tp) 



(16) 



4.1 Solvability of Problem (16) 

Since (—00, +00) x (0, T] is unbounded, we first confine our attention to the truncated 
version of (16) in a finite domain Qn = (—N,N) x (0, T]. Let u N (x,9) be the solution (if 
it exists) to the following problem 



' min {d e u N - Cu N , u N - (e x - K)+) = 0, (x, 9) G Q N , 

' N , 9€(0,T\, 

x G (—N,N). 



d x u N (-N,9) = 0, d x u N (N,9) 
^ u N (x,0) = (e x -K)+, 



(17) 



In order to prove the existence of solution to problem (17), we construct a penalty 
approximation of problem (17). Suppose (x,9) satisfies 



f d e u? - £u? + P £ {u^ - 7r £ (e x - K)) = 0, (x, 9) G Q 



AT, 



d x u»(-N,9) = 0, d x u?(N,9) 



where [3 e (t) (see Fig. 1.) satisfies 

Pe(t) € C 2 (-00,+00), P E {t) <0, 

^(t)>0, %(t)<0, &(0) = -C , 
where Co > is to be determined. Note that 

'0, t > 0, 



9 G (0, T], 
x€(-N,N), 



(18) 



lim/3 e (t) 

£->0 



and that vr £ (t) (see Fig. 2.) satisfies 



-00, t < 0, 
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and 7r e (t) G C°°, < ir' e (t) < 1, <'(t) > 0, lim n £ (t) = t A 



T 




-Co 



Fig. 1. p e {t) Fig - 2 - 

Lemma 4 There exists a unique solution u^(x,9) E Wp' (Qjv) to problem (18) for any 
p > 1. Moreover, the following estimates hold, 

ir £ (e x -K)< v%{x, 9)<k9 + x 2 + e x+ ^- pc ^ +e + 1, (19) 

< ^uf (x, 9) < e x+ ( b -P c V +e , (20) 

d e u?(x,9)>0, (21) 



where k > max < c + 



2 , (b-pcA-ic 2 ) 2 n 



27(l-p 2 )c : 



^ ' _i_ 
; C -t- 



7 (l-p 2 )c 2 e( i >-" cA ) +T -(6-pcA-ic 2 )] : 



2 7 (l-p 2 )c 2 



2 1 

Proof: The Leray-Schauder fixed point theorem implies the existence of W p ' solution to 
problem (18), and the comparison principle holds for the strong solution. The proof of the 
uniqueness is standard. 

First, we prove estimate (19). Note that u\{x,9) := 7r e (e x — K) satisfies d$ui — Cu\ + 
/3 e (ui — 7r e (-)) < in Qn, if we choose 



1 



(22) 



ui(x,0) 



-&(0) = C = pc\e» + - 7 (1 - ^)cV JV - 

Moreover, when N is large enough, we have d x u\{— N, 9) = 0, d x u\{N,9) 
uf (x, 0). Thus u\{x,9) = ir £ (e x — K) is a subsolution to problem (18). 

On the other hand, one can check u^ix, 9) := k9 + x 2 + e x +( b ~P cX ) 9 + 1 satisfies super- 
solution property dgU2 — Cu 2 + f3 £ (u2 — ^e(')) > in Qn. Moreover, we have 

d x u 2 (-N, 9) = -2N + e -N+(b-pc\)+e < Qj 

^n 2 (iV, 9) = 2N + e ^+(0-p C A)+e > g iv 

u 2 (s,0) = x 2 + e x + l > 7r E (e x — K). 

Applying the comparison principle, we conclude (19). 

Next, we verify inequality (20). If we differentiate the equation in (18) w.r.t. x, then 
v(x, 9) = d x u^ (x, 9) satisfies 



d e v - \c 2 d xx v - (b - pcX - \c 2 )d x i 
+ 7 (1 - p 2 )c 2 vd x v + P' E (-){v - Ti> £ e x ) = 0, (a?, 9) G Q 

v(-N,9) = 0, v{N,9) 

{ v(x,0)=7r' £ (-)e x , 



IN, 



-,N 



9e(0,T], 
x e (—N,N). 



(23) 
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Since v\ = and V2(x,9) = e x+ ( b P cX ) +d are subsolution and supersolution of (23), the 
estimate (20) follows by comparison principle. 

Denote u s (x, 9) := u*(x, 9 + 5) for < 5 < T, u s (x, 9) satisfies 

' d e u 5 - Cu 5 + p e (u s - K e (e x - K)) = 0, (x, 6) G (-N, N) x (0, T - 5], 

< d x u s {-N,9) = 0, d x u s {N,6) = e N , 9 G (0,T - 6], 

s u 5 {x,0) =u*(x,S) >7r £ {e x -K) = u^{x,0), x€(-N,N). 

Combining with (18), applying the comparison principle, we have 

u 5 (x, 9) > uf (x, 9), (x, 9) G (-N, N) x (0, T - 5], 

which implies (21). □ 

Lemma 5 Problem (17) has a unique solution u N (x,9) G Wpj oc (Qiv) ^C(Qn) satisfying 

(e x - K) + < u N {x, 9)<k9 + x 2 + e ^+(b- P cX)+0 + i, (24) 
< d x u N (x, 9) < e x+(b-pc\)+e^ (25) 
d e u N {x,9)>0, (26) 

where k is defined in Lemma 4. 

PROOF: Let Cn be a generic constant independent of e. 

Since {x, 9) > ir £ {e x — K), then \/3 £ (u^ — ir £ (e x — K))\ lp iq n \ < Cn- One can treat the 
nonlinear term in the equation in (18) as a linear term with the coefficient ^7(1 — p 2 )c 2 d x u^ . 
Thanks to (20), applying C a,a / 2 estimate (De Giorgi-Nash-Moser estimate, [13]) with a = 
1/2 , we have 

\ u e lc* 1 /2,i/4( Qjv ) < C N - 

Letting e — > 0, we have a continuous limit (of a subsequence if necessary) u N (x,9) up to 
the boundary, that is, 

u?(x,9) ^u N (x,9) in C(Q N ). 

2 1 

In view of (20), applying W p ' estimate [9], we have 

\ u ?\w?\Q*r\Bs) ~ ° N ' 

where B$ is a disk with center (In K, 0) and radius 5 > 0. Hence u N (x, 9) G W^ oc (Q.zv) and 

(x, 9) - u N (x, 9) weakly in W*'JjQ N ), 
which also implies u N (x,9) is a 

W P,Lc solution of (17). Furthermore, (24)-(26) are conse- 
quences of (19)- (21). 
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2 1 

In this below, we show the uniqueness. If not, we can find two W p ' solutions u\ and 112 
satisfying (24) and (25), and Af = {(x,9) G Qn '■ u\ > U2} ^ 0. Observe that, dgu\ = Cu\ 
and dgU2 > Cu2 in Af. Thus, u\ — U2 satisfies 



\c 2 )d x {ui - u 2 ) 



dg(ui - u 2 ) - \c 2 d xx {ui - u 2 ) - (b - pcX 2 , 
+i 7 (l - p 2 )c 2 {d xUl + d x u 2 )d x ( Ul - u 2 ) < 0, {x, 9) G Af, 

(ux-u 2 ){x,e) = 0, (x,6) G d p Af\{x = ±N}, 

d x (u 1 -u 2 )(x,9) = 0, (x, 9) G dpAf f) {x = ±N}. 

Owing to (25), applying maximum principle (see [17]), we have u\ — U2 < on Af, which is 
a contradiction to the definition of Af. □ 

Lemma 6 There exists a unique solution u{x,6) G W£ loc (Q) n C(Q) to problem (16) sat- 
isfying 

(e x - K) + < u(x, 6)<k9 + x 2 + ^+{b- P c\)+e + ^ ( 2 7) 

< d x u(x, 9) < (?+(b-pc\)+e^ ( 28 ) 

d e u(x,9)>0, (29) 

where k is defined in Lemma 4. 

Proof: First, we claim problem (17) is equivalent to the following problem 

,N _ (e* _ #)} = 0, (x,6)EQ N , 



mm{dQU N — Cu ±r 



N 1 



-N,6) = 0, 8 x u N (N,9) 



^ u N {x,0) = (e x -K)+, 



9e(0,T], 
x G (~N,N). 



(30) 



In fact, suppose w{x,9) is a W ' loc solution to problem (30), the maximum principle [17] 
implies w > 0, combine with w > e x — K to get w(x,9) > [e x — K) + . Hence, if is also a 
solution of (17). Together with the uniqueness of (17), the equivalence follows. 
Now we can rewrite problem (17) as 



' d e u N -Cu N = f(x,0), 
8 x u N (-N,9) = 0, 8 x u N (N,9) 



,v 



(31) 



(x,9) G Q N , 
e", 0<E(O,T], 
u «(x,0) = (e x -K)+, xe(-N,N), 

where f(x, 9) = X{u N {x,e)={e x -K)+} ( ~~ (b — pc\)e x + ^7(1 — p 2 )c 2 e 2x ) . Thanks to (25), if we 
apply Wp 2,1 interior estimate [9] to (31) for arbitrary M < N, 

\ uN \w£\Q M \B s )^ C M, 

where B$ is a disk with center (In K, 0) and radius 5. We emphasize Cm only depends on 
M, but not on N. Fix M > 0, and let TV — > +00, u N leads to a limit (possibly a 

subsequence) in the fixed domain Qm in the sense, 



u 



N 



U 



(M) 



weakly in W*) (Qm) as N ->■ 00. 



p,loc 
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Moreover the sobolev embedding theorem implies 

U N U {M) in c{Q M ), d x u N -> d x u^ i n C (Qm). 

It is clear that u(x, 9) := u^ M \x, 9), (x, 6) £ Qm is well defined and u is the solution to 
problem (16). Moreover, d x u £ C(Q) and we can deduce u £ C(Q) from the C a estimate. 
(27)-(29) are consequences of (24)- (26). Lemma 1 and 2 implies the uniqueness of solution 
to problem (16). □ 

4.2 The Obstacle of (16) 

Now we will characterize the optimal exercising time of V. Problem (16) is an optimal 
stopping problem, which gives rise to a free boundary that can be expressed as a single- 
valued function of 9. For later use, we define 

S ■= {(x, 9) £Q: u(x, 9) = ie x - K) + ] (Stopping region), 
C := {(x, 9) £ Q: u(x, 9) > (e x - K) + } (Continuation region), 

where u is the solution of (16). Thanks to the continuity of u from Lemma 6, S is closed 
and C is open, respectively. 

Lemma 7 There exists S(x) : (— oo, +oo) — > [0, T\ such that 

S = {(x, 9) £ Q : < 9 < S(x)}. (32) 
Proof: If (x ,0 ) £ S, i.e., 

u(x ,9 Q ) = (e X0 -K)+, 

according to (27) and (29), we have 

u(x o ,0) = (e xo - K) + , 9£(0,9 }. 

Hence {xo} x (0,9o] C 5, so we can define S(x) : (— cxo,+oo) — > [0, T] as 

[ 0, if u(x,6) > {e x -K) + for any 9 £ (0,T], 
\ sup{9 £ (0, T] : u(x, 9) = (e x - K)+). 

By the definition of S(x), (32) holds. □ 

Lemma 8 The free boundary S(x) is strictly increasing w.r.t. x on {x : < S(x) < T}. 

Proof: As in the proof of the equivalence between (17) and (30) in Lemma 6, we can show 
that problem (16) is equivalent to the following problem 

J mm {d e u-Cu, u-(e x -K)} = 0, {x,9) £ Q, 
{ u(x,0) = (e x - K)+, x£R. 



S(x) 
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When u(x, 9) = (e x — K) + = (e x — K), we have 



; - K > 0, 

)u(x,9)-£u(x,9) > 0, 



which implies 



^-H*$r^j?>- (34) 



For any xq G {x : < S*(x) < T}, denote 5(xo) = #o £ (0, ^Z 1 ) , in view of (32), we know 
u(x , 9) = e x ° — K > 0, € (0, e }. 
Denote Qo = (— oo, +oo) x (0, #o] an d define a new function u(x, 9) on Qo by 
_ / ^ j u(x,0), (x,9) £ (-oo,x ] x [0,9 ], 
" 1 e*-if, (z,0) G [x ,+oo) x [0,e ]. 



U[X, 



Since {xq} x (0,9q] <5) then u(x,9), d x ll(x,9) are continuous in Qo- Now we want to 
prove u(x, 9) is the solution to problem (16) in the domain Qo- 
By the definition of u{x, 9) we know 

u(x,0) = {e x -K) + , 

According to (34), we can check u(x,9) satisfies 

min {d e u-Cu, u- {e x - K) + } = 0, (x,t) G Q - 

Hence, we know u(x,9) is the solution of (16) in domain Qo- By the uniqueness of 
W P,ioc(Q) n C (Q) solution which satisfies (27)-(28) to problem (16) we know 

u(x,9) =u(x,0), (x,0)eQq. 

In particular, 

u(x, 9) = u(x, 9) = {e x - K) + , (x, 9) G [x , +oo) x (0, 9 ]. 

By the definition of S(x) we know 

S(x) >9 Q = S{x q ), x > Xq, 

therefore the monotonicity of S{x) is proved. 

Next we will prove the strict monotonicity of S{x) on the set {x : < S(x) < T}. 
Suppose not (see Fig. 3). There exists x\ < X2 such that 

S{x) = 9 G (0,T), xe{x u x 2 ). 

Thus 

f d e u(x,9) -Cu(x,0) = 0, (x,0)e(x 1 ,x 2 )x(0 o ,T], 
1 u(x,0 o ) = (e x - K), xe(x 1 ,x 2 ). 
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Hence 



d e u\ 



1 



0=0o 



7 (1 - p 2 )c 2 e 2x + {b- pc\y 



< 



1 



(b - pcX) - ^7(1 - p 2 )c 2 e x 



1 



< -glC 1 " P )c'e^ + (6 - pcA)e Xl < 0, 
which contradicts (29), the two inequalities in (35) is due to (34). □ 



(35) 




Xl x 2 

Fig. 3. No strict monotonicity of S{x) 




XQ 



■ X 



Fig. 4. Discontinuity of S(x) 
Lemma 9 The free boundary S(x) is continuous on the set {x : < S(x) < T}. 
Proof: Suppose not (see Fig. 4). There would exist an xq such that 

S{x Q ) := 9 1 > hm_S(x) := 6 2 - 

Then we would have 

d e u(x Q , 6) = d x9 u(x , 9) = 0, Be (9 2 , 9 X ). 
Since dgu > and in the domain (xq — e, xq) x (62, 0\) 

de(d e u) - ^c 2 d xx (d e u) - [b - pcX - X -c 2 - 7(1 - p 2 )c 2 d x u\ d x {d e u) = 0. 

Applying Hopf's Lemma [5] we obtain 

d x8 u(x , 9) < or dgu = 0, (x, 9) G (x - e, x ) x (9 2 , 0i), 

which results in a contradiction. □ 

Since S(x) is strictly increasing with respect to x on the set {x : < S(x) < T}, 
then there exists an inverse function of S(x), we denote it as s(9) = S~ 1 (x), < 9 < T. 
Note that the strictly monotonicity of S(x) is equivalent to the continuity of s(9), and the 
continuity of S(x) is equivalent to the strictly monotonicity of s(9). Owing to Lemma 8 
and Lemma 9, we give the following theorem. 

Lemma 10 There exists an optimal exercising boundary s(9) : (0,T] — > R such that 

S = {(x,9) € Q: x> s{9)}. (36) 
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Moreover, s(9) is strictly increasing with respect to 9 and 



s{9) > z , (37) 
a(O) := lim s(9) = x , (38) 

( InK, b- pc\< 0, 

where xq = < 2(b-ocX) ^ n particular, dgu is continuous 

y max{mK, In ^_pi^2 j> b — pcX > 0. 

across s{9) and s(9) € C[0, T] n C°°(0, T]. 
Proof: We can define 

5 _1 (a:), < 9 < T, 

inf{x : S(x) =T}, 9 = T. 



8(6) 



Equation (36) is the consequence of the definition and monotonicity of S(x). According to 
(34) in the proof of Lemma 8, we know (37) is true. Lemma 8 implies s(9) is monotonic 

increasing with respect to 9, then we can define s(0) := lim s(9), and the proof of (38) is 

0— >o+ 

similar to the proof of strictly monotonicity of S(x) in Lemma 8. 

Moreover, since dgu > and (e x — K) + is a lower obstacle, from [6] we know dgu is 
continuous across s(9) and s(9) € C°°(0, T\. □ 

4.3 Main Result: Characterization 

Now, we are ready to present the complete characterization of the value function V and the 
indifference price P. 

Theorem 11 1. The indifference price of (7) is independent of the initial capital w, 
and P(w,y,t) := P(y,t) is the unique W^f oc (Q 2 ) flC(Q 2 ) solution of the variational 
inequality (13)-(14) satisfying \d y P(y,t)\ < C for some constant C. 

2. The value function of (6) has the form of V(w,y,t) = Ut(w + P(y,t)), and is the 
unique solution of (9)-(10) in W^f ^.(Q}) (~lC(Q 1 ) satisfying 

\d y V(w,y,t)\<C\V(w,y,t)\. (39) 

3. There exists C°° function y* : [0, T) R, such that the strategy defined by 

< = -ez.Y v d y P{Y v ,v) + ±, ue(t,r% 

t* = inf{z/ > t : Y v > y*(u)} A T, 
is optimal of the control problem (6), in the sense that 

V(w,y,t)=E[U 7 *(WZ + g(Y^ y ))\T t ], 
where W t * = w + J t 1 -K^aiXdv + dB v ). 
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Proof: By the fact yd y f(y,t) = 8 x u(x,t) and estimates of (28), we conclude \d y f(y,t)\ < 
C for some constant C. Applying Lemma 6, Lemma 2, together with Lemma 1 in order, 
we conclude (l)-(2) of Theorem 11. Thanks to the representations of (12) and (15), the 
optimal control can be written as 

f < = -£ ■Y v d y P{Y v ,v) + ±, v € (i, t*), 
\ r* = mf{u > t : (W*,Y V , v) £ C[V}} A T. 

Note that the continuation region of (8) satisfies 

C[V] = {(w, y, t) € Q 1 : V(w, y, t) > U t (w + g(y))} 

= {(w, y, t) € Q 1 : U t {w + P(y, t)) > U t (w + g(y))} 
= {(w,y,t) G Q} : P(y,t)) > g(y)}. 

Therefore, the obstacles of P and V are the same, and the optimal stopping time r* can be 
written invariant to W*: 

t* = hh> > t : P(Y U , v) < g(Y u )} A T. 

Take y*(t) = e s ^ T ~ l \ where s(-) is the function in Lemma 10. Thanks to the the result of 
Lemma 10, together with the transformation between P(y,t) and u(x,t), we conclude the 
representation of r*. □ 



5 Dual Representation and Ramifications 

In this part, we will give an alternative proof of the result on dual representation of the 
indifference price given in Proposition 7 of [10]. Thanks to the characterization Theorem 11, 
the proof based on the stochastic control approach is rather straightforward. 

Now we recall the market with stock prices given by (1) and (2). Let Z be the collection 
of all the processes of the form 

Zf = exp | - ^X 2 t - A£?i j exp | - ^ J ip 2 v dv - j Lp u dB u ^, 

where <p is some T v progressively measurable process satisfying (p^dv < oo. Then, the 
collection of equivalent local martingale measures A can be written by 

A = {Q^ : dQ^ = Z*dF, eZ}. 

Among of them, we refer Q° to the minimal martingale measure (MMM). Under Q^, two 
processes 

B$ := B t + Xt and Bf := B t + [ ip u dv (40) 

Jo 

are both standard Brownian motions. Let Qf and Ft be the probability measures of 
and P restricted to Tt- Then, we have 

dQf = ZfdF t . 
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To proceed, we introduce the following concept. The relative entropy of a probability 
measure Pi with respect to P2 is defined by 



ff(Pi|P 3 ) 



OO, 



log 



Pi « P 2 , 
otherwise. 



Example 1 The relative entropy with respect to P is 



£T(Q*'|P) = E ( 

= E' 
1 



log- 



tfP 



1 



-A 2 r - \b t 



(pldu 



ip u dB v 



AT + E 8 - 



Example 2 Given a stopping time r £ 7o,T, the relative entropy Qr with respect to is 



H( 



E^ 



E^ 



E^ 



log- 
1 

~ 2 



E^ 



log 



Lp v dB 



ifldv 



Proposition 12 TTie indifference price P(w, y, t) of (7) admits following dual representa- 
tion for all (w, y) elxK + , 



P(w,y,0) = esssu PreroT essinfQ^ A 1x9? \g(X?> )] + -H(Qf\Q° T 



(41) 



Proof: Thanks to Example 2, it is enough to show that, for all t £ [0,T] 



P(w, y, t) = J(y, t) := esssup^T^essinfpgs <j E^ g(Y^ ) + — / ^efc/ 



where $ is the collection of all T v progressively measurable process satisfying l p 2 dv < 00, 



and En.]=E<H- l-^l- 



Recall that in (40), both B x and 5^ are Brownian motions under Q v . We can rewrite 
the process Y of (1) in terms of B x and B^, 



dY u 



cp\ - cy/1 - p 2 (p u )Y u dv + c P Y v dBl + cy/1 - p 2 Y v dBt 
cpX - ca/1 — p 2 ip v )Y v dv + cY u dBv ,ip , 



(42) 



where is another Q^-Brownian motion. 
Since Q 1 ^ ~ P for arbitrary <p> £ <J>, 
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Jq ip 2 dv < oo P-almost surely 

is equivalent to 

Jq <p 2 du < oo Q^-almost surely 

and vice versa. In other words, there exists an 1-1 map r : i— > <I>, such that the distribution 
of T((p) under P is equal to the distribution of ip under Q^. Define 



dY? = {b- c P X - cy 7 ! - f?<p v )Y*du + cY^dB v . (43) 

Comparing (43) with (42), we conclude Y under is equal to Y r ^ under P in distribution. 
Therefore, we write J(y, t) as a standard control problem of the following form, 



J(y, t) = esssup Te7 ; T essinf^ U t \g{Y^^) + i- jf r^)^] | 
= esssup reTt T essin% e $ |lE t g{Y^ y ' t ) + ^- tp 2 v dv .1 



(44) 



The second equality of (44) follows from the fact = r($). 

Applying exactly the same procedure of the verification theorem Lemma 1, we can 
conclude J(y,t) = v(y,t), provided that there exists v G W^f oc (Q 2 ) solves 

min I - d t v - \c 2 y 2 d yy v - (b - cp\)yd y v + sup^ \cy/l - p 2 y<pd y v - ^V? 2 }, 

v(y,t)-g(y)}=0, (y,t)eQ 2 , ^) 

[ v(y,T)=g(y), y G R+. 

Thanks to Theorem 11(1), P(w,y,t) = P(y,t) G W^ oc (Q 2 ) solves (13)-(14). By utilizing 
the quadratic structure of sup^j • } in (45), one can check P(y, t) also solves (45) . Therefore, 
P(y,t) = v{y,t) = J(y,t). □ 

How does the price P(y, t) change, if we scale its payoff g by n times, or if we change 
the risk aversion parameter 7? As a result of the dual representation Proposition 12, we 
have the following properties on the price P(y,t). 

Proposition 13 P(y, t) decreases with respect to 7 and A, and increases with respect to b, 
satisfying 

nP[g](y,t) > P[ng](y,t) (n > 1), (46) 

where P[<p\ stands for the indifference price with the payoff function ip. 

Proof: Note that, Proposition 12 remains valid for the P[ng], if we change g into ng in 
(41). Since H(Qr\Qr) ^ s non-negative, the monotonicity in 7 and the non-linearity of (46) 
follow directly from Proposition 12. 

Suppose Yi (i = 1,2) are two processes of (43) related to Aj). If b\ > b<i and Ai < A2, 
then Y\ > Yi almost surely by comparison result of stochastic differential equation, and 
(44) implies Pi > P2, where Pi are the associated prices. □ 
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An alternative proof of Proposition 13 using a PDE approach is given in the Appendix, 
and it is interesting as its own right. In fact, Proposition 13 reveals natural economic 
facts. For instance, since 7 can be interpreted as the absolute risk aversion coefficient, 
proposition 13 claims that the employee's risk preference directly affect his exercise behavior. 
It implies a more prudent agent (with a larger coefficient of risk aversion 7) would be likely 
to exercise the option earlier to realize a cash benefit, and then invest it in other asset to 
earn the time value of the money obtained from the exercise of the option. So when he 
exercises, he gets less value than the risky agent, hence the value function decreases w.r.t. 
7. In addition, (46) implies that if the agent owns n (n > 1) pieces of options, compare 
with owning one piece of the option, the agent care less about the value of every piece of 
the option. So the average value of the indifference price about these pieces of options is 
less than the indifference price about one piece of this option. 

As a straightforward consequence of (27) and (28), we present the following estimates 
on P(y,t) and V(w,y,t). 

Proposition 14 P(y, t) and V(w, y, t) satisfy 

{y - K) + < P(y,t) < k(T -t) + (my) 2 + ye (b- P cX)+(T-t) + ^ (47) 
0<d y P(y,t)< e {b ~ pcX)+ (T - 4) , (48) 

_ e -l(w+(y-K)+)+±\ 2 t < V < _ e -7(^+fc(r-*)+(lBy) 2 +2/e( 6 -^) + ( T -*)+l)+|A 2 t ) (49) 

< d y V(w, y, t) < - 7 e( b -" cA ) + ( T -'VK y, t), (50) 
where k is defined in Lemma 4- 

6 Application to ESO Costs 

In this part, we briefly describe the application of the above result to Employee stock options 
(ESO). An ESO is a call option on the common stock of a company, issued as a form of 
non-cash compensations. Compared to the American call, the main differences of ESO are 
the vesting period and job termination risk. 

Suppose the company's stock price follows Y of (1), which is non-tradable in the market. 
Given a unit of ESO with maturity T, strike K, and vesting period t v € (0, T), its payoff 
is equivalent to the conventional American call only if the exercise time occurs between t v 
and T, otherwise zero. Therefore, we can write its payoff as 

(Y T - K) + I {T > tv} 

for the exercise time r G To t- I n addition, if the job termination is taken into account, we 
have the following revised payoff, 

(Y rAr a — K) + I{ TAr a> tv }, 

where r a is the time of employee's job termination. We assume r Q is a random variable 
having exponential distribution with parameter a > 0, which is independent of Brownian 
motions B, and B. 
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As suggested by FASB rules, the price of ESO is evaluated by risk-neutral measure, 
under which stock price Y is a martingale. For simplicity, we assume b = in (1), and P 
is the risk- neutral measure specified above. Following the arguments of [11, Section 5], the 
indifference price has a representation given as follows. 

1. When ESO survives throughout the vesting period, that is, if t > t v , then the ESO 
cost C(-) satisfies 

C(y,t) =EQ[(Y T , Ara -K)+] 

= E<2[e- a ( T * "'Hi;* - K)+ + J7* a{Y u - K)+ e - a ^du], 

where r* = inf{^ > t : Y v > y*(u)} A T is optimal exercise time in Theorem 11. This 
corresponds to the following PDE characterization: C(y, t) is the unique C 2,1 solution 
of 

8 t C + \c 2 y 2 d yy C - aC + a{y - K)+ = 0, (y,t) G Q 2 n {y < y*(t)} n {t v < t < T}, 
C(y, t) = (y- K)+, (y, t) G Q 2 n {y > y*(t)} n {t v < t < T}, 
C(y,T) = (y-K)+, y E M + , 

(51) 

since the boundary curve y*(t) is smooth due to Theorem 11 and PDE is non- 
degenerate locally. 

2. If t < t v and ESO is alive at the moment, then the ESO cost is 

C(y,t) =WP[C{Y tv ,t v )I {Ta>tv} ] =EQ[e~ a ^C(Y tv ,t v )}. 
Thus, C(y,t) follows, 

d t C+^c 2 y 2 d yy C-aC = 0, (y,t) € M + X [0,t v ), 

with terminal condition C(y,t v ) given by the solution of (51). 

Acknowledgement We thank Professor Nicolai V. Krylov for the discussion on the gen- 
eralized Ito's formula. 

A Appendix 

A.l A Generalized Ito Formula 

We present a generalized Ito formula to the functions in Sobolev spaces which is an ex- 
tension of the result in [8]. We thank Prof. Krylov for confirming this result by email 
communication. 

Proposition 15 Let X be a diffusion on the filtered probability space (Q, J-, P,F = {Ft}), 
with generator L, initial time s, and initial state x € Suppose v 6 i oc (Q) for some 

open set Q C Define 

T s ' x = inf{r > s : X r <£ Q}. 
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Then, for any ¥ -stopping time r < Tq S , we have 

E[v(X t ,t)} = v{x,s) +E / Lv(X r ,r)dr 



(52) 



(53) 



Proof: For any F-stopping time r < Tq X , by Theorem 2.10.2 of [8], we have 

E[v(X T , t)] < v(x, s) + E I Lv(X r , r)dr 

L is 

If we apply (53) to a function v\ = — v, then 

E[v 1 (X t ,t)]<v 1 (x,s)+E J L Vl {X r ,r)dr 
which yields 

E[v(X t ,t)\ > v(x,s) +E / Lv(X r ,r)dr . (54) 

J s 

From (53) and (54), we conclude equality holds for (53). □ 
A. 2 Proof of Proposition 13 

One can conclude the result of Proposition 13 by the following two lemmas. 

Lemma 16 The solution u(x, 0) to the problem (16) decreases with respect to 7 and A, and 
increases with respect to b. 

Proof: We may prove all three monotonicities in the same way by comparison principle. 
Hence we only present the proof of the monotonicity of u(x, 9) with respect to 7. Suppose 
7i > 72) and u^^(x,9)(i = 1, 2) is the solution to the following problem 



- hc 2 d xx u N 



pc\ - \c 2 )d x u^ 



+i 7 i(l - p 2 )c 2 (d x u^) 2 + /3 £ (u^ - 7r e (0) = 0, (x, 9) e Q N 
d x u^(-N,9) = 0, d x u? {i) (N,9) = e N , 
u^ ) (x,0)=Tr £ (e*-K), 



Set w(x,9) := u^Jx, 0) — u3 2 Jx,6). Then w(x,9) satisfies 



ee(o,T], 

x € (-N, N), 



dew - \-c 2 d xx w - (b - pcX - ]rc 2 )d x w + -71(1 - p 2 )c 2 (d x u^ + d x u^ 2) )d x w 



1 



+#(■)«» = 2^ " Ti)(l " pVO^)) 2 < 0. 
Combining with the initial and boundary conditions, we have 

u^ 1) (x,9)-u^ 2) (x,9)=w(x,9) <0. 
Letting e — > 0, N — > +00, we know u(x, 9) is decreasing w.r.t. 7. □ 
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Lemma 17 The solution to problem (16) satisfies 

nu[g(x)] > u[rig(x)] (n > 1), (55) 

where g(x) = (e x — K) + , and u\g\ represents the solution to problem (16) with the obstacle 
and initial condition g. 

Proof: Set u(x,9) := nu[g(x)], then u(x,9) satisfies 

minion — Cu — \^n{n — 1)(1 — p 2 )c 2 (d x u[g]) 2 , u — ng} = 0, 
u(x, 0) = ng(x). 

Note that u[ng(x)] satisfies 

min{(dgu — Cu)(x, 9), u — ng} = 0, 
u(x, 0) = rig(x), 

We can confine the above problems in the bounded domain Qn- Suppose u^f^] and u N [rig\ 
are the solutions of the following problems 

' min {d e u N [g\ - Cu N [g\ - ^n(n - 1)(1 - p 2 )c 2 (d x u N [g\) 2 , 

u N [g\-ng(x)} = 0, (x,9) 6 Q N , 

* d x u N [g\(-N,9) = 0, d x u N [g\(N,9) = ne N , 9 € (0,T], (56) 

k u N [g\(x,0)=ng(x), X e(-N,N). 

min {dgu N [rig\ — Cu N [rig\, u N [rig\ — rig(x)} = 0, (x,9) € Qn 
< d x u N [ng\(-N,9) = 0, 8 x u N [ng\(N , 9) = ne N , G (0,T], (57) 
u N [ng](x,0) = ng(x), x € (-N,N). 

Comparing (56) with (57), letting N — > oo, we obtain (55). □ 
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